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An  interdisciplinary  approach  blending  systems  theory,  numerical 
analysis,  and  mathematical  software  was  employed  to  address  various  aspects  of 
the  numerical  solution  of  matrix  Riccati  equations  and  some  related  problems. 
A  variety  of  numerical  issues  related  to  the  computational  solution  of  the 
algebraic  matrix  Riccati  equation  were  studied.  Our  approach  concentrated  on 
the  generalized  eigenproblem  formulation  for  the  solution  of  general  forms  of 
algebraic  Riccati  equations  arising  in  both  continuous-  and  discrete-time 
applications.  These  general  forms  result  from  control  and  filtering  problems 
for  systems  in  generalized  (or  implicit  or  descriptor)  state  space  form. 

The  issue  of  numerical  condition  of  the  Riccati  problem  was  addressed  as 
was  the  use  of  various  types  of  balancing  — both  numerical  and  system- 
theoretic —  to  improve  numerical  condition.  A  rather  sophisticated  Fortran 
software  package  called  RICPACK  was  developed  for  use  both  as  a  research  tool 
and  as  a  production  vehicle.  Numerical  experiments  with  RICPACK  were  per¬ 
formed  to  investigate  a  number  of  proposed  condition  numbers.  None  of  these 
numbers,  both  those  developed  under  this  project  and  those  developed  by 
others,  have  been  found  to  yield  completely  satisfactory  results.  The  subject 
remains  open  still.  However,  overall  experience  with  RICPACK  to  date  indi¬ 
cates  that  it  is  the  most  powerful  software  yet  developed  to  solve  general 
classes  of  Riccati  equations  reliably  for  modest-sized  (say,  of  order  a  few 
hundred  or  less)  problems. 


Particular  attention  was  focused  in  the  latter  stages  of  the  project  on 
the  special  structure  of  models  of  physical  systems  represented  in  linear 
second-order  form.  Exploiting  that  structure  in  solving  associated  Riccati 
equations  was  studied.  As  a  hy-product,  convenient  numerical  tests  for 
controllability  and  observability  were  derived  directly  in  terms  of  the 
original  second-order-model  matrices. 
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